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Abstract: We study linear cosmological perturbations in a previously introduced family 
of deformations of general relativity characterized by the absence of new degrees of freedom. 
The homogeneous and isotropic background in this class of theories is unmodified and is 
described by the usual Friedmann equations. The theory of cosmological perturbations is 
modified and the relevant deformation parameter has the dimension of length. Gravita- 
tional perturbations of the scalar type can be described by a certain relativistic potential 
related to the matter perturbations just as in general relativity. A system of differential 
equations describing the evolution of this potential and of the stress-energy density per- 
turbations is obtained. We find that the evolution of scalar perturbations proceeds with 
a modified effective time-dependent speed of sound, which, contrary to the case of general 
relativity, does not vanish even at the matter-dominated stage. In a broad range of values 
of the length parameter controlling the deformation, a specific transition from the regime of 
modified gravity to the regime of general relativity in the evolution of scalar perturbations 
takes place during the radiation domination. In this case, the resulting power spectrum of 
perturbations in radiation and dark matter is suppressed on the comoving spatial scales 
that enter the Hubble radius before this transition. We estimate the bounds on the de- 
formation parameter for which this suppression does not lead to observable consequences. 
Evolution of scalar perturbations at the inflationary stage is modified but very slightly and 
the primordial spectrum generated during inflation is not noticeably different from the one 
obtained in general relativity. 
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1. Introduction 

The future astrophysical and cosmological data will, among other things, allow one to 
perform deeper tests of our theory of gravity, general relativity (GR). Although there are 
no compelling reasons to question seriously the validity of GR, several puzzles associated 
with its long-distance behavior, primarily those connected with the phenomena of dark 
energy and dark matter, make it an important task to see whether GR continues to be as 
good a description of gravity on astrophysical and cosmological scales as it is on the scales 
of Solar System. Modern literature is abundant with various proposals of modified gravity 
addressing this issue. A common feature of almost all such modifications is that new 
propagating degrees of freedom are introduced in the gravitational sector, in the form of 
additional fields or higher derivatives of the metric field. A typical example is f(R) gravity, 
which in its original variables is a purely metric theory with higher derivatives but also can 
be transformed into a scalar-tensor theory. Extensions of this kind typically change both 
the Friedmann equations describing the evolution of a homogeneous and isotropic universe 
and the equations describing cosmological perturbations, by which modifications they can 
be tested. 

It is, perhaps, less known that one can modify GR without introducing any new prop- 
agating degrees of freedom. Modification of this kind, in particular, should be expected 
not to change the homogeneous background evolution in cosmology, affecting only the way 
the perturbations evolve. Some of the generic features of this scheme were explored in Q 
without specifying a concrete theory. 

A specific class of modifications of GR without new propagating degrees of freedom 
was proposed by one of us in 0], with the same class of theories envisaged earlier in a 
different formulation in ||. These theories legitimately can be called deformations of GR 
in the sense that there exist continuous paths in the parameter space smoothly connecting 
any of them with GR. The theories in this class provide a new framework for testing general 
relativity, since, using the observational data, one can determine how close is the real-world 
gravity to GR within a class of its deformations. In one of our previous papers [|], we were 
concerned with the spherically symmetric solution in this class of theories, its features and 
possible implications. In the present paper, we turn to the description of cosmology. Our 
aim is to set a stage for possible tests of the theories of |2| by working out the modified 
equations describing the evolution of cosmological perturbations. 

Before describing the class of theories under investigation in detail, it is worth dis- 
cussing their basic properties in elementary terms. If viewed as purely metric theories, 
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they can be regarded as the usual effective theories with a nonlocal Lagrangian that can 
be expanded into an infinite series of local terms built from the curvature invariants: 
C = aiR + 02-R 2 + a^RfivR^ + a^R iiVaT R ilv<yT + ... The coefficients in front of these 
invariants, however, are not arbitrary in the present case, but are such that the whole infi- 
nite series can be resummed with the help of auxiliary non-propagating fields. In the end, 
one obtains a local second-derivative theory with only two propagating degrees of freedom 
(two polarizations of the gravitational wave), and a certain limit can be taken to recover 
GR if one so wishes, see ||. 

To elucidate this idea, we illustrate it on a simple example of a single field with a 
higher-derivative action. Let us start from the linearized Lagrangian 

£i 2 )[h] = hiUh, (1.1) 

where □ = V^V^ and our signature convention is (— , +, +, +). One can think of the field 
h as that of the graviton with the factor of the inverse Newton's constant G^ 1 absorbed 
into the definition of h to endow it with the canonical mass dimension [h] = 1. All tensor 
indices of h have been discarded in order to make the discussion as simple as possible. 

Lagrangian ( |1.1[ ) can be modified by introducing higher-derivative non-renormalizable 
terms. For any finite number of such terms, the resulting theory has additional propagating 
modes corresponding to the poles of the arising propagator. However, when an infinite 
number of such higher-derivative terms is present, they sometimes can be resummed into 
a denominator so that the only pole of the arising propagator is still that of the massless 
graviton. One of the simplest examples realizing this idea is 

c[h\ = \h{u + t 2 u 2 + e 4 n 3 + £ 6 n 4 + ...)h = \h l ^ l2u h ■ (1.2) 

Here, £ is the length scale of the modification. It is clear that the usual "graviton" is the 
only propagating mode in (|1.2|) . To see this explicitly, note that the non-local theory (|1.2| ) 
can be obtained from the local one 

C[h, X } = ^x{0-m 2 )x + xOh+hnh, (1.3) 
where m = £~ l , by integrating out the auxiliary field x- Indeed, the equation for x that 



stems from flL3J) is formally solved by 



Dh , , 

m z — Li 



Substituting this into (1.3), one obtains ( |1.2| ). At the level of action ( |1.3[ ), the limit m — > oo 



that gives back (LI) works by making the field x infinitely massive, and thus effectively 
decoupling it from the theory. At the same time, writing action (|1.3| ) in terms of the new 
field variable h = h + x, 

C[h, X } = hah-^m 2 x 2 , (1.5) 



one can see that the field x is an auxiliary non-propagating field, and the theory (|l.3| ) has 
the same number of degrees of freedom and structure as the one we started from, Eq. (|1 . 1|) . 



- 3 - 



Thus, the described modification ( |l.2| ) is just a field redefinition, which gives another way 
to understand why no new propagating degrees of freedom is introduced in this case. 1 

The described elementary modification scheme is (almost) identical to the one at play 
in the class of theories Q at the linearized level. However, in theories it is extended to 
a full non-linear level in a non-trivial fashion. It then turns out that, by the same trick of 
introducing auxiliary non-propagating fields, quite non-trivial terms in the arising effective 
metric Lagrangian can be reproduced, for instance, the Riemann 3 term of importance in 
quantum gravity at two loops; see ||. 

To summarize, the idea of the deformations of GR is that higher-derivative terms in 
the effective metric Lagrangian can be added without introducing extra degrees of freedom, 
and the theory is then made local at the expense of introducing some non-propagating 
auxiliary fields. In the class of theories this idea is realized elegantly by combining the 
metric and auxiliary fields into a single field — an SU(2) Lie-algebra-valued two-form. 

Another general comment should be made about the class of theories under investiga- 
tion in this paper. The above discussion seemingly implies that they can be regarded as 
ultraviolet modifications of general relativity. Indeed, usually, modifications generated by 
the addition of higher-derivative terms to the Hilbert-Einstein Lagrangian are of relevance 
at high energies, e.g., close to the Planck scale, but are unimportant at low energies of rel- 
evance in cosmology. However, our above example shows that this interpretation may not 
always be valid. Indeed, the theory fll.2j ) has the field and dynamical content equivalent to 
that of (1.1). This property will be lost after truncation of the series in (1.2) or after the 



violation of the delicate relations between the coefficients in this infinite series. Through 
these relations, the high-derivative terms (potentially important in the ultraviolet) are all 
connected with the low-derivative terms (important in the infrared). In the case of the 
modified gravity theory to be considered, this UV/IR interplay will be seen in the fact that 
the theory can be generated from the self-dual formulation of general relativity simply by 
making the cosmological constant a function of the variable which in GR has the meaning 
of the Weyl curvature (more on this below). This allows gravity to be modified on cosmo- 
logical scales rather naturally, simply by making the usual cosmological constant a slowly 
varying function of the curvature. At the same time, it is easy to remain consistent with 
GR on Solar-System scales by requiring the "cosmological function" to be approximately 
constant in the relevant region of curvatures (see the next section for more details). In 
this sense, the theory simultaneously looks like an infrared modification of gravity. This 
interplay between the ultraviolet and infrared modification is an interesting feature that 
characterizes the class of theories under investigation. 

The plan of the paper is as follows. In Sec. ||, we briefly describe the family of defor- 
mations of GR which is the subject of this study. In Sec. ||, we describe the homogeneous 
isotropic Universe and the usual Friedmann equations in the language of two-forms adopted 
in this paper. Section || describes the linearized field equations around the evolving cos- 



ine cannot discard the field x i n (1-5) if the "matter" part of the action interacts with the original 
field h = h — x- In this case, there arises a theory with non-propagating auxiliary scalar x that couples to 
matter, resembling the modified source gravity Q. This will also be the case in the modified-gravity theory 
under consideration in this paper. 
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mological background. In Sec. || we classify the types of perturbations and describe the 
physical gauge-invariant quantities. Section ^ obtains equations describing perturbations of 
the scalar type. These are analyzed in details in Sec. ^. We study the effect of modification 
on the CMB and matter power spectrum in Sec. ||. We conclude with a brief discus- 
sion of the results obtained. Preliminary results concerning tensor modes are described in 
Appendix |B|. 

2. A family of deformations of GR 

In this section, we briefly review the theory under investigation in this work. For more 
details, the reader is directed to a recent description in 0. 

2.1 Preliminaries: two-form field and the metric 

As we mentioned in the introduction, in our theory the metric and the auxiliary non- 
propagating fields are "unified" into a single two-form field B^ u , i = 1,2,3, where n,v are 
spacetime indices. The (complexified) rotation group SO(3, C) acts on the two-form 
via B^ v — > OjBftv, Oj £ SO(3, C), and the theory is invariant under the corresponding 
gauge transformations as well as under the action of spacetime diffeomorphisms. 

The physical metric of the spacetime, to which all fields are supposed to couple uni- 
versally, is determined uniquely by the two-form field. To begin with, we require that 
the triple of two- forms B l be self-dual with respect to this metric, which determines its 
conformal class. In other words, introducing the associated volume form e^p^ for a metric 
g^y and imposing the self-duality conditions 2 

\^ pa B; a = iB^ (2.1) 

fixes the metric g^ u modulo a conformal rescaling. 

An equivalent explicit description of the conformal class of metrics determined by B l ^ u 
is given by the formula 

yj-detig^g^ oc e a ^ 5 B^B^B^ , (2.2) 

where e a/3/yS is a densitized completely anti-symmetric tensor having components ±1 in any 
coordinate system, and the proportionality symbol is used here to denote that the metric is 
defined up to conformal rescalings. The reality conditions imposed on the theory are that 
the conformal metric ( |2.2D is real Lorentzian. To fix a unique metric from its conformal 
class, it suffices to specify its volume form. Below, we shall explain how the triple B 1 ^ 
defines not only a conformal class of metrics, but a volume form as well. 

Vice versa, given a real metric g^ u with Lorentzian signature, one can introduce the 
corresponding tetrad one-forms 91, I = 0, 1, 2, 3, so that the line element reads 

ds 2 = 9 I ®9 J r]u, (2.3) 
2 The spacetime indices in (|2.l[) are raised with the use of the inverse metric ' . 
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where rjjj is the Minkowski metric. Then, introducing an arbitrary time plus space split 
I = (0, a) of the "internal" index /, we can define a set of self-dual metric two-forms T, a : 

S a = i6° A 9 a - ^ e abc 9 b A 9 C . (2.4) 

Here, the letters a, b, c are new SO(3) indices. It is convenient to distinguish between the 
original SO(3) indices in B % and the new index that appears in the metric two-forms S a , 
which justifies the notation. Let us also give an expression for a convenient basis in the 
space of anti-self-dual two- forms: 

S a = i6° A6 a + ^ e abc 9 b A 9 C . (2.5) 

Any other set of self-dual two- forms can be decomposed in the basis provided by S a . 
Thus, we can write 

B i = b^ a . (2.6) 

Therefore, the main dynamical field of the theory can always be chosen as a collection of a 
metric g^ that, in turn, defines the metric forms S a , together with the nine scalars forming 
the matrix b\. These scalars will later be seen to be non-dynamical and determined by 
other fields present in the system. We do not impose any independent reality conditions 
on b\ as such reality conditions are induced once b\ are determined in terms of other 
fields. Thus, b % a are in general complex fields. Note that, as the metric g^ u undergoes 
a conformal transformation g^ v — > Q 2 g^ I/ , the metric self-dual two- forms transform as 
S a —7- f2 2 £ a . Thus, if the scalars b % a transform as b l a — > Q,~ 2 b % a: then the original two-forms 
B l are unchanged. A certain normalization condition on the scalars b l a will be introduced 
below to fix this conformal freedom, thus fixing a particular metric from the conformal 
class of (^). 

2.2 The gravitational action 

The vacuum theory of gravity under consideration is described by the following action: 



S[B,A] 



8vrG 



B l AF\A)--V{B i AB J ) 



(2.7) 



Here, A i is an SO(3,C) connection field, F^A) = dA i + (l/2)e ijk A^ A A k is its curvature, 
and y(M y ) is a (holomorphic) function of a complex symmetric 3x3 matrix variable M lJ . 
It is required to be a homogeneous function of degree one: V (aM y ) = aV (M Jjf ) for any 
number a, so that it can be applied to a matrix-valued B % A B 3 four-form, with the result 
being again a four-form; see 0] for more details. The function V(M^) is also required 



to be a scalar with respect to the SO(3, C) gauge transformations. In action (2J), this 
function plays the role of potential for the B l field. 



Using representation (2.6), let us introduce the "internal metric" 
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which is a symmetric 3x3 matrix, in general complex. It turns out to be useful to 
parameterize it by its trace h = h l3 5ij and trace- free part H 13 : 

h ij = K (gij + . (2.9) 

Then the potential V (B i A B 3 ) gives rise to the SO(3, C)-invariant function V (h 13 ) , which 
can be presented as follows: 

V (hv) = ±h U(H), U(H) = A - ^2 TrF 2 + O {H 3 ) . (2.10) 

Here, U(H) is a function of the traceless matrix H 13 with dimensions of curvature, and 
we have expanded it in powers of H %3 assuming that this potential function is analytic in 
the neighborhood of H %3 = 0. Below we will see that the case of h %3 = S 13 corresponds to 
general relativity, so that the traceless matrix H 13 parameterizes the deformation away from 
GR. The constant term Ao in decomposition ( |2.10| ) will play the role of the cosmological 
constant. The quantity I is a new parameter of the theory with dimensions of length. 
In principle, in the theory under consideration, there is scope for making this parameter 
complex, but, in this paper, we shall only consider the simplest case of real I. The numerical 
constant in the expansion in ( 2.1C| ) is introduced for future convenience. We will see that 



only negative sign of the second term in ( 2.10| ) describes a theory without instabilities, 
which explains our choice of this sign. 

2.3 Coupling to matter 



As we have stated in Sec. |2T , the physical metric which universally couples to matter 
belongs to the conformal class determined by the two- form field B % . It is necessary to fix 
the remaining conformal freedom and to specify this metric uniquely. 

Under conformal transformations, described in the end of Sec. |2.1| , the "internal met- 
ric" introduced in ( |2.8[ ) transforms as h %3 — > O One can specify the physical metric 
by fixing the conformal ambiguity for h %3 . To do this, we introduce another scalar (holo- 
morphic) function R (h 13 ) of the matrix h ij that is required to be homogeneous of degree 
one in the components of h %3 : R (a/i iJ ) = aR (/i Jjf ). The conformal freedom of h %3 is then 
fixed by the condition 

R (h ij ) = 1 . (2.11) 

While this prescription for fixing the physical metric may seem arbitrary at first, it can be 
shown to arise quite naturally by considering the motion of a "small body" in the theory 
described by (2/7); see @. 



We normalize the function R (M 3 ) so that R (8 13 ) = 1. The potential R (h i3 ) can then 
be decomposed similarly to ( 2. 10| ) : 



R(h) = h U m {H) , U m (H) = 1 - |Tr# 2 + O {H 3 ) , (2.12) 

where U m (H) is a new function of the traceless matrix H = H 13 , again assumed to be 
analytic in the neighborhood of H 13 = 0. The quantity g is a dimensionless parameter. 
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Just as with the case of £, in the theory under consideration, there is scope for making the 
parameter g complex. However, in this paper, we shall only consider the simplest case of 
real I and g. Note that, once the condition ( |2.11| ) is imposed, one can use it to express the 
trace part of the matrix h %3 in terms of its traceless part. 

We shall see that only the terms written in ( 2.1C| ) and ( p. 12 ) are going to matter for 



the linearized theory we are going to consider. Thus, the only new parameters that we will 
have to consider on top of those available in GR are the length parameter I in ( |2.10| ) and 



the dimensionless constant g in ( 2.12 ) that can be of any sign. 

In the class of theories under consideration, matter fields will in general couple to 
the two-form B % that plays the role of the fundamental geometrical field of the theory. 
However, for the purposes of this paper, we assume that our matter content is "standard" 
in the sense that our material components couple to the spacetime metric in the conformal 
class defined by B % and further selected by the condition fl2.11|) . Thus, we assume that the 
dependence of the matter action S m on B % is such that S m [B] is a functional of g^ u and 
depends on the scalars b % a solely via the function R(h). This allows us to write: 

5S m = J ^-det(^) [r^gT + ^f|^a) • ( 2 -!3) 

However, since the two-form field B l does not transform under the conformal transfor- 
mations bg^y = eg^ v , 5b' l a = — any material action that arises from S m [B] is also 
conformally-invariant. This immediately gives: 

OSrndRi 

1 ~ dR db\ a ' [ > 

where T := T^ u g^ v , or, using the homogeneity of R(h): 

dS m T 



dR 2R 

In other words, the variation of the material action contains the following two terms: 

T dR 



(2.15) 



5S m = I A /-det ( 5m ,) ( T^SgT + ^r-q^K ) . (2.16) 



'a 



After the condition (|2.11| ) is imposed, the quantity T^y becomes the usual stress-energy 
tensor of matter. We can anticipate already at this stage that the field equations will 
contain terms involving the first derivative dR/db l a of the potential R(h). 

2.4 Field equations 

The set of equations obtained by varying the full action of the theory (given by ([2]7]) plus 
the material action) with respect to A 1 (assuming that the material action does not depend 
on A 1 ) is 

dB i + e m^(B) A B k = . (2.17) 

These equations are solved by 

A ^ {B) := J^B^^^^ ' d6t B := -Y/ 3kB ^ BlPBk p" ■ (2 - 18) 
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We note that this connection is conformally-invariant and has the correct transformation 
properties of an SO(3) connection. That is, when the two-form field B* transforms as 
8B^ V = e l i k u)j B^ vl the connection transforms as S^A 1 ^ = u J A? l — d li uj i . A demonstration 
of this involves a simple identity satisfied by the two-forms B^ u ; see ||]. We also note that, 
in most cases, instead of using ( 2,1§| ) for finding A l (B), it is easier to solve ( [2.17 ) directly. 



Once the connection A l {B) is found, one computes its curvature: 



F i := dA i + ^ e ijk A j A A k . (2.19) 
The field equations then read, in the parametrization fl2.6|): 

b\F [A(bH)] = (Jjtffti + ^Abibfj S b + 2vrGT {^b% + -A^&f) S b -2vrGT ab S fc . 

(2.20) 

Here, A and A m are the Legendre transforms 

BU BTI 
A = U{H) -dm Hl3 > ^m = U m (H)-^-HV (2 . 21) 

of the potentials U(H) and U m (H), respectively, T is the trace of the stress-energy tensor 
of matter, T ab are the components of the traceless part of the stress-energy tensor, and S a 
are the anti-self-dual metric two- forms ( |2.5[) . For the matter in the form of ideal fluid, we 
have 

T=(p-3p), T ab = {p + p)[ Sa f + 2ie abc = = )- 3S ab , (2.22) 

where p and p are the energy density and pressure, respectively, u a is (proportional to) the 
momentum vector, and S ab is the traceless matrix describing shear. 

As anticipated in the previous subsection, the right-hand side of our field equations 
includes a term proportional to the trace T of the stress-energy tensor and containing the 
first derivative of the potential function R{h) (or U m (H)) of the material sector. The 
remainder of the dependence of the right-hand side of ( |2.20| ) on the stress-energy tensor 
components can be easily understood by working out how matter couples to gravity in the 
so-called Plebahski formulation of GR; see, e.g., ||. The only change in the present case as 
compared to GR in the Plebahski formulation is in the appearance of the derivative terms 
dU m (H)/dH . A detailed derivation of the field equations ( |2.20 ) can be found in M. 



In ( |2.21 ) , we have introduced the "cosmological function" A as the Legendre transform 



of U(H). When H 13 = 0, the above field equations reduce to those of GR in Plebahski 
formulation. For later purposes, we note that the function A is, in fact, a function A(^) 
of the quantity 

dU 

$ = $y := (2.23) 

that can be related to the usual Weyl curvature (at least for small deviations from GR). 
It is equally possible to parameterize the theory either by the original potential U{H) or 
by its Legendre transform A(\I>). However, as we shall discuss in more detail a little below, 
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the nature of the modification becomes more clear in the A(^) parametrization. This is 
due to the fact that the field H lJ that measures the departure from GR arises in this 
parametrization as the first derivative 

= --^L (2.24) 

of the cosmological function A. Thus, for a function A(\I / ) which is constant in its domain 
to a high degree, the deviations from GR are small. 

2.5 Large-distance modifications 

As we have already mentioned in the introduction, the class of theories under investiga- 
tion can be viewed as ultraviolet modifications of gravity. This is most clearly seen in 
the parametrization of the theory by the potential U(H). When it is expanded around 
H l i = 0, the first non-trivial term contains a mass parameter m 2 = l/£ 2 , and modifica- 
tions away from GR only manifest themselves at energies E larger than m. The low-energy 
limit E <C m of any of our theories is given by GR. However, when the theory is looked 
at non-perturbatively, e.g., in the Hamiltonian formulation (see ||), one finds that the 
parametrization by the Legendre transform A(\I/) of U(H) is more natural. More specifi- 
cally, one finds that the departure of the theory from GR is connected with the fact that 
the cosmological constant becomes a function A(^) of what in GR used to be the Weyl 
part of the Riemann curvature tensor. The parametrization by A(\&) then allows one to 
arrange for modification of gravity at any desirable scale of curvatures while keeping it 
close to GR at the scales of the Solar System. Indeed, the theory behaves close to GR 
whenever the first derivative (2.24) of the function A('I') is small. Thus, to be consistent 



with the Solar-System tests of gravity, one should assume the function A(\&) to be approx- 
imately constant in the range of curvatures probed in the solar neighborhood. One can 
then allow A( X I / ) to vary in the range of curvatures much smaller than those in the Solar 
System, which will produce deviations from GR at large distances (see Q as well as a more 
accessible exposition flCl ). Thus, according to our modification scenario, the value of the 
cosmological constant in the Solar System can be different from the value Ao measured in 
cosmology, with the function A(\P) interpolating between these two values. The parameter 
I arising m (ETIcD is then given by the second derivative of the function A(\l/) at *$> tJ = 0. 
Since gravity is not tested directly at very large distances, there are only very weak direct 
constraints on this parameter. As we shall see, a study of the physics of the cosmic mi- 
crowave background (CMB) and formation of the large-scale structure as predicted by this 
theory will place much stronger constraints on the value of t. 

Another important point that we would like to make is in the following. A detailed 
analysis shows that, as the first derivative ( |2.24| ) of A(^) becomes of order unity (in the 
sense that any of the eigenvalues of the symmetric traceless matrix H l;) become of order 
one), the theory under consideration develops a specific degeneracy, which results in a 
singularity of metric d2.2|) Q]. It can also be argued that it is no longer consistent to use 
the classical theory in this regime because quantum effects become important. Thus, it is 
necessary to avoid this regime if one wants to remain in the domain of applicability of the 
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classical physics, which is usually the case for large distances of relevance for cosmology. 
Therefore, the classical theory will remain applicable only when the condition 



OA 



< 1 (2.25) 



is satisfied. This condition can be rephrased by saying that the departure from GR as 
measured by the quantity ( 2.24| ) must be small. We assume this condition in what follows. 



2.6 The cosmological constant and vacuum energy 

In general relativity, the cosmological constant A can be regarded as corresponding to 
the energy density of the vacuum, and one is free to absorb it in the energy density and 
pressure or into the scalar-field potential V(<p) in all equations without exception. That is, 
the transformation 

A -)■ A — A, + P^P-^t; (2-26) 



8vrG" r c 8vrG 



in the case of ideal fluid or 



A^A-A, V&) + — (2.27) 

in the case of scalar field leaves invariant all equations of the theory. 

In the theory under consideration, this property is true only for the equations describing 
the homogeneous and isotropic universe, which are the same as in general relativity. It is not 
true in the general case, which can be seen in ( |2.20| ), in which the trace of the stress-energy 
tensor of matter T enters differently from the genuine cosmological constant Ao [present in 
the "cosmological function" A(\P)]. This feature of the theory of modified gravity persists 
in the equations for cosmological perturbations, as we shall see below. An exception is the 
special case when the parameter £ 2 is fixed to be £ 2 = 1/4<?Aq. In this case, the invariance 
with respect to ( 2.26j ) is restored at the level of linear perturbations. 



3. Homogeneous and isotropic universe 

Consider a homogeneous and isotropic background described by the two- forms Bq. On 
such a background, the traceless matrix H lJ describing the departure from GR is equal to 
zero due to the spacetime symmetries. Thus, we have Hq = 0, and we can choose the 
scalars b l 0a = 8 l a , so that Bq = Sq. The background equations are then the Friedmann 
equations in the Plebahski formalism: 

pi [A(Eo)] = 27rCTo + A E , _ 27rG - q ^ (31) 

Here, 

T :=(p-3p), f :=(p + p) (3.2) 

are the components of the background stress-energy tensor. We emphasize that the back- 
ground evolves in our theory in exactly the same way as in GR. Modifications only arise 
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in the dynamics of perturbations. This is, however, the story for an ideally homogeneous 
backgrounds. In reality, there are inhomogeneities on various scales that need to be aver- 
aged over. It is natural to expect that this will introduce modifications of the background 
equations as well, to be taken into account. However, since the averaging issue is not fully 
understood even in the simpler case of GR, we postpone such an analysis to our future 
publications, and, in this paper, study the cosmological perturbations around an ideally 
homogeneous and isotropic background. 

In the conformal time rj, the metric of a spatially flat homogenous and isotropic universe 

is: 



dsQ = a 2 (rj) 



(3.3) 



In what follows, we often omit the argument of the scale factor a{rf). 

A set of self-dual and anti-self-dual two-forms describing this metric can be chosen as 

4 = a 2 (id V A dx i - ^e ijk dx j A dxA , ^ = a 2 (idr] A dx i + ^e ijk dx j A daA . (3.4) 

The connection A l compatible with the set of self-dual two-forms £q, i.e., satisfying dT, l + 
e ljk Al A Eq = 0, is given by 

Al = mdx i , (3.5) 

where 

%:=-, (3.6) 
a 

and the prime denotes the derivative with respect to the conformal time rj. The curvature 
F l = dA l + \^ k A^ A A k of connection fljT5| ) is given by 



where we have used the relations 

d V A dx i = (4 + £ l o) , dx i A dx j = ^ ijk (sjj - E§) (3.8) 

to decompose Fq into its self-dual and anti-self-dual parts. 

Substituting ( |3.7D into ( |3.1[) , we obtain the dynamical equations for the background. 
The condition that the self-dual and anti-self-dual parts of the field equations are separately 
satisfied gives two equations whose symbolic linear combinations £ — £ and 3S + £ are 
just the Friedmann equations 

U 2 8ttG P A 2n' + U 2 o ^ . , on , 

^ = -T + T' — = -&rGp + Ao. (3.9) 

4. Linearized field equations 



Linearizing the field equations ( [2.20 ) around the homogeneous and isotropic background, 



we denote the perturbation of a quantity by a symbol 5 next to the symbol representing 
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it. There arises a perturbation of the metric two- forms 5T, 1 as well as a perturbation of 
the scalars Sb l a . Using the fact that Bq u = 8 l a , now we can identify the two types of SO (3) 
indices. Thus, from now on we drop the distinction between the two types, and use only 
the i,j,k indices. The condition R{h 1 ^ = 1 implies, at the first order in perturbation, 
that the matrix 5b lJ is traceless. The perturbations of the metric h l;) and of its trace- free 
part are then 5W = SIP* = 25b&\ Since Hq = 0, the first-order perturbation of the 
quantity TrH 2 that appears in the potential functions is zero. Thus, 5A = and 5A m = 0. 
After some algebra, the linearized field equations are found to be 

S^+Sb^F^ = —6TZi-2irG5T ij f:i+ 27TGTo + A ° > g-27rGT ,5g-47rG/afr ij E J . (4.1) 
3 3 

Only the last terms on the left-hand side and on the right-hand side of this equation are 
not present in GR, with all other terms being exactly like in the Plebahski formulation of 
GR. The perturbation of the components of the stress-energy tensor of fluid is given by 

ST = 5 P - 35p , 5T 1 ' = (Sp + 5p)5 ij + 2\e ijk 5u k - 35S ij , (4.2) 

where, for convenience, we have absorbed the factor p + p into the momentum vector u\ 
The quantity 5F l is the linearized curvature corresponding to the perturbation of the B l 
field given by 

8& = 8b is Y? Q + 8& , (4.3) 
where 5Y? is the perturbation of the metric self-dual two-forms (see below for explicit 



expressions). In (4.1), we have introduced a new parameter k with dimensions of the 
energy density: 

k := (g - I) T + - = fi{p - 3p) + 1PA , (4.4) 



3 J 8ttG \£ 2 3 

where To = p — 3p is the trace of the material stress-energy tensor, p\ = Aq/8ttG is the 
energy density of the cosmological constant, and i and g are defined in (|2~10| ) and (gig ), 



respectively. The quantities (3 and 7 are two new convenient dimensionless parameters: 

These parameters have the following meaning. The parameter 7 determines the deviation 
of the gravitational potential V (h 1 -*) from one given by V (h 1 ^) = Ao-R This case 

corresponds to 4/3 = 7 and describes the theory that couples in the same way to the 
cosmological constant and to the trace of the stress-energy tensor. The parameter (5 controls 
the deviation of the physical metric from the so-called Urbantke metric which coincides 
with ( |2.2[ ) precisely. The value /3 = corresponds to coupling of matter to the Urbantke 
metric. Finally, the limit to GR is obtained by keeping j3 fixed (and arbitrary), and sending 
I — > (thus 7 — > 00) keeping the product ^5B % ^ fixed. To summarize, schematically we 
have: 

(3 = <^=> coupling to the Urbantke metric , 7—7-00 <?=^ limit to GR . (4.6) 
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Now, using the expression ([T?]) for the background curvature and equations ([O]), we 

get 

2ttGT + A H 2 + H' ,= H 2 -H' 

3 = ' 27rGT ° = "ST" • (47) 

The linear equations for perturbations can then be written in the form 

SB ij Fr> + 6F* - A(^ 2 + n')6& + ^{U 2 - H')8& 
2a z 2a z 

= -^-<5TE - 2irG5T 3 S J - AttGk 5B l] YP Q , (4.8) 

where the stress-energy perturbations of the fluid are given, as before, by 

ST = [5p - 38p), ST 1 ' = (s p + Sp )S^ + 2\e ijk 5u k - 35S ij , (4.9) 

and 5B lJ is the non-metric part of the perturbation (see the next section for explicit expres- 
sions). This form of the linearized equations is most convenient for further calculations. 

5. Classification of perturbations 



A general infinitesimal perturbation of the background (3.4) can be decomposed into the 
background two- forms: 

SB 1 = 6$ ij & + 5T ij f% . (5.1) 

However, it is clear that much of this perturbation is a gauge transformation. Indeed, 
the anti-symmetric part describes an infinitesimal SO (3) rotation of the triple of 

two-forms E l and is thus a pure gauge. Moreover, some part of the general perturbation 



(5.1) is an infinitesimal diffeomorphism of the background (3^4). To separate the physical 
quantities from the gauge ones, we thus need to compute the effect of an infinitesimal 
diffeomorphism on the background fl3.4| ). The Lie derivative of B i , suitably corrected by 
a gauge transformation, can be shown to be given by C^B 1 = T>B l (£), where B l (£) is the 
interior multiplication of a vector field £ with two- form B l , and T> is the covariant derivative 



with respect to the -B-compatible connection. For the background (3.4), a straightforward 
but lengthy calculation gives 

= ~ [mf + f + ac] s + £ M £ ' + + G) + & + O)' + s fe 

+ \ ie - ac] e* + + c,ij] ^ + \<> ik m + oy - n . (5.2) 

Here, we have decomposed the vector field £ into its components £ = + (A, where 

£i is transverse f j = 0. We see that the effect of a diffeomorphism is a change in the 
metric described by X together with a gauge transformation on E . The latter can always 
be corrected by an SO(3) rotation, which allows us to drop the anti-symmetric part of 
the tensor <53>* J from now on. Note that diffeomorphisms do not affect the symmetric 
traceless part of 5<$> l i , i.e., they do not act on the non-metric 3 components of the B i field. 



3 We shall often refer to the non-dynamical scalars present in B l in addition to the metric as its "non- 
metric" components. 
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In principle, this could have been expected, for they act in the spacetime manifold, not in 
the internal space, where the non-metric components of B % reside. 

Let us now consider the usual classification of perturbations into scalar, vector and 
tensor sectors and deduce the most general form of these perturbations after all gauge 
freedom is fixed. 



5.1 Scalar sector 

Scalar perturbations are described by scalar functions, so that we have 

SQ* = + 4j [x,ij ~ + , (5-3) 

ST*i = #^ + 1 - ^5ij A<j\ + ^ k oj >k . (5.4) 

Here, we have introduced the factors 1/a 2 for future convenience, (j> and ip are the trace 
parts, x an d o are the scalar parts of the traceless symmetric parts, and and oj are the 
scalar parts of the anti-symmetric parts of the self-dual and anti-self-dual perturbations, 
respectively. It is easy to see that diffeomorphism ( |5.2| ) can be used to set the scalars a 
and u) to zero, while an SO(3) gauge transformation can be used to set to zero the scalar 
6. Thus, we are led to consider the following "physical" scalar perturbations: 

5B i = SB^l + + ^ , (5.5) 

where we have introduced a convenient notation for the non-metric part 

5B ' 13 = ^2 U,ij ~ l^jAx) . (5.6) 

In other words, the scalar perturbation ( |5,5| ) is a sum of a non-metric perturbation SB^'Eq 
and the usual metric perturbation 

5^ = + if>f? . (5.7) 



The potential x can ! i n general, be complex, while (p an d ip are real. The potential x 
does not have a counterpart in GR, while the two potentials <p an d ip have a very simple 
relation to the standard GR potentials appearing in the scalar perturbation of the metric: 



ds 2 = a 2 



(1 + 2$ GR ) drj 2 - (1 - 2* GR ) (dx'Y 



(5.8) 



An elementary calculation shows the relation between $gr> ^gr an d the potentials (f), tp 
in (O): 



$GR 



3-0 + <j> 



GR 



1p-<p 



(5.9) 
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5.2 Vector sector 



The most general perturbation of the vector type can be reduced to 



6B % 



C^Sj + ^Sjl , (5.10) 



where the vector potentials C and C are transverse, = C,i = 0) an d C ma y be complex. 
Alternatively, it is a sum of a non-metric perturbation 5B l3 Y? Q with 

5BV = 2C (i ' j) (5.11) 

and a metric one: 

££* = 2£to)Ej . (5.12) 

5.3 Tensor sector 

Perturbations of the tensor type that correspond to gravitons are given by 

6B i = ^X ij ^ j + p ij V j , (5.13) 

where both x*- 7 and p l i are symmetric, traceless and transverse (x^ ,i = P l \i = 0) matrices, 

and x %3 ma Y be complex while p %3 is real (see, however, Appendix). Again, the perturbation 
is a sum of a non-metric one 

SB* = l x « (5.14) 

cr 

and a metric one 

(JE^p^'Sj. (5.15) 



6. Scalar perturbations 
6.1 Curvature 

In this section, we consider the important case of scalar perturbations in detail. First 
of all, we compute the connections and then (linearized) curvatures for the perturbations 
described by 4>, ip, and \. To compute the perturbation of the incompatible connection for 
a perturbation SB 1 , one has to solve the algebraic equations Vq(5B % ) + e ljk SA J A Sq = 0, 
where £q is the background two-form. The curvature of an infinitesimal connection 5 A 1 is 
then given by F' l (5A) = Vq5A 1 , where Vq is the covariant derivative with respect to the 
background connection. The computations are straightforward, but rather lengthy. The 
expressions for the connection components are: 



8A\ 
5A 



SA{ 



3i 



ip,idr) 



(aV)' 



Sij - e 



1 

2i I 

-^Ax,idr] - -e ijk Ax,kdx J 



ijk(P,k 



dx J 



{ a *^ - ^ A * 6 ij) dxJ 



(6.1) 
(6.2) 
(6.3) 
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The expressions for the curvatures, decomposed into their self-dual and anti-self-dual 
parts, are 



1 



+ 



V5AI 



1 
1 

1 

4^2 



-2^- + <%(</>" + A<j>) + 2ie 



ijk 



4>,k 

a 



a — 2H(f>'5i 



E> , (6.4) 



+ 2iej,-& 



4^,y + <5y« - A^) + 4% f «V + + 
2^,o- + SijW + Atp) + 4<% (u'lj) - 2H 2 ip + 



^0 



(6.5) 



1 

2a 4 " 
1 

2^4 



(d 2 - 2Wfy - A 



8 H A 



4^ 2 ) 
" 3% Ax 



<%A X 



(6.6) 

where (51?* J in the last formula is given by (5.6), and Fq is the background curvature (|3.7|). 
The differential operators in brackets act only on \ but not on the background forms Sq. 

6.2 Field equations 

Now we can substitute the expressions for the curvatures obtained above into the linearized 
field equations ([4.81). For the scalar sector, the sources are 



5T = 5p- 35p, 5T ij = {5p + 6p)5 ij + ^e ijk 5u k - 3 \6S tij - ^ASs) , (6.7) 

where we have introduced a convenient normalization of the scalar part of the velocity 
perturbation 6u. The resulting equations are most usefully separated into the symmetric 
traceless, antisymmetric, and trace parts of the matrices arising in front of the background 
two- forms £q and Eg. Let us first analyze the symmetric trace- free parts. After removing 
the directional derivatives and multiplying by a 2 in the self-dual sector, and by 2a 2 in the 
anti-self-dual sector, the self-dual and anti-self-dual parts give, respectively, 

1 



^ + 2a 2 k" " mX ' ~ AX ~ = 4nGKX ' 

<p + ip + 



1 



1 



X 



-Ax 



(6.8) 
(6.9) 



-12irGa 2 5S. 

a" \ 3 

It is more convenient, however, to consider these equations as those determining the usual 
relativistic potentials ^gr and ^gr; see (|5.9|). We have 



*GR := 
$GR 



*P-4> = 
2 

(j> + 3ip 



1 (I 



A-kGkx + -s -Ax + Ux + 2^ 2 x + QTrGa 2 5S , 



(6.10) 



AttGkx + ^ (-x" + \^x + n X ' + 2n 2 X ^j - 67rGa 2 5S . (6.11) 
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The difference between the relativistic potentials is given by 

* GR - $GR = ^ U" ~ + l2TrGa 2 5S . (6.12) 

Thus, even in the absence of shear, the relativistic potentials ^gr and $gr are not equal 
to each other in our theory, which is typical of modified gravity. 

The antisymmetric part is only present in the anti-self-dual sector. After removing the 
directional derivative and multiplying by ia 2 , one obtains 

*' GR + H$ GR - ^ A X ' = AirGaSu , (6.13) 

which, after solving for the perturbations ^gRj ^GR) and X) wu l determine the scalar part 
of the 4-velocity perturbation. Finally, the trace parts of the arising equations are most 
convenient in their schematic linear combinations S — S and 3S + 2. Using the equations 
of motion for the background, after some simple algebra, we get 

A* G R - 3W GR - 3H 2 $ G R - ^2 A 'x = ^Ga 2 5p , (6.14) 
* GR + 3W GR + {W! + U 2 )^ G R + ~A($ GR - *Gr) + ^($GR " *Gr)' 

+ {2W + U 2 ){$ G R - ^gr) - ^2 A 'x = ^Ga 2 5p . (6.15) 

Thus, the arising equations are exactly like in GR apart from the presence of simple addi- 
tional terms containing x- 

6.3 Conservation equations 

The following standard conservation equations hold in the theory under consideration: 

Sp' + m{8p + Sp) = 3{p + p)y' GR + ^A5u , (6.16) 
5u + 3H5u = a5p + a{p + p)<S> GR + aASS . (6.17) 

This can be verified directly using equations ( |6.9D , ( |6.10| ), and ( |6 . 1 3| ) above. Note that the 
field X; responsible for the modification, does not appear in these equations. This is of no 
surprise as the stress-energy is conserved in this theory in the usual way. 

6.4 Equations for the "non-metricity" scalar x 

At this point, the usual GR procedure would be to use the fact that ^gr = ^GR (in 
the absence of shear) and obtain a system of equations for one of the potentials. In our 
case, this is not possible because any such equation will involve contributions containing 
the nonmetric scalar x- However, it is possible to obtain a closed system of equations for 
this scalar. Thus, we solve for <3?gr and ^gr i n terms of the new potential x and deal 
with differential equations for x- The resulting equations are not simple, but it is worth 
presenting them here to show that their structure is similar to the corresponding equations 
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in GR. Thus, we now substitute ip + (f> = $gr — ^GRj ^gr> and ^GR from ( &9j) , (|6.1C| ), 
and ( |6.11 ), respectively, into ( |6,14| ) and ( 3.15 ), and, after some lengthy algebra, obtain 



Q-TJ 

h 2 a x - —x (W + n 2 ) + m 2 in 2 - m') x 



4irGa 2 5p = A(AttGk X ) - 3U {4ttGkx)' - 3U 2 {4ttGk X ) 
2 

- H- \\ - ■ 

2 

+ 6nGa 2 [ASS - 3U (5S' + HSS)] , 
4TiGa 2 5p = {AttGkx)" + 3U {A-kGkx)' + {W! + V 2 ) (^Gk X ) 



(6.18) 



+ 



n 2 x" + \a x (h' - u 2 ) + \x (n" + m'n - 3H 3 ) 



+ x (2H"n + 2%' 2 - 2n 2 n' - n 



5S" + hU {SS' + U5S) - -ASS 



(6.19) 



+ 6nGa 2 

Similarly, substituting ^gr and ^GR into ( |6.13| ), we get 

AirGaSu = (4ttGk X )' + H (^Gk X ) + -j [x (W + + ^ {2U 1 - U 2 )] 

+ 6vrGa 2 (SS' + HSS) . (6.20) 

This equation determines the scalar part of the velocity perturbation once the perturbations 
X and SS have been found. 

Now, introducing the speed of sound c s via 



Sp = c s Sp + t5£ , 



(6.21) 



we can write a single dynamical equation for x with the entropy perturbation S£ as the 
source. It reads: 

^Ga 2 T5£ = (AttGkx)" - c 2 s A {AttGkx) + 3ft(l + c 2 s ) (^Gk X )' 
+ [2W + U 2 (l + 3c 2 )] (AttGkx) 

+ ^(h 2 x" + ^x [n'-u 2 (l + 3c 2 )] 



+ y [%" + H'H (7 + 3c 2 ) + 3H 3 (c 2 - l)] 
+ x [2U"U + 2U' 2 - 2U 2 W (1 - 3c 2 ) - U A (l + 3c 2 )] J 
SS" + (5 + 3c 2 ) H (SS' + HSS) - Q + c 2 J ASS 



6vrGa 2 



(6.22) 



In the absence of shear {SS = 0), this dynamical equation can be solved for x> after which 
one can compute all other quantities of interest, such as the relativistic potentials and 
perturbations of energy density and pressure. 
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6.5 New relativistic potential 

The obtained complicated equations for x can be considerably simplified by introducing a 
new quantity 

$ := UkGk + ^\ X ■ (6.23) 

One rationale for considering precisely this combination is the identity that follows directly 
from |[I|) and ( |6720| ) : 

5p + —5u - -A5S J = A$ . (6.24) 

Precisely the same identity holds in the case of GR with on the right-hand side of this 
equation being the relativistic potential $GR) which suggests that the quantity <I> should 
play the role of the main relativistic potential in our theory. Below we shall see that this 
expectation is realized. 

Anticipating the role that is going to play, it is illuminating to rewrite the main 
dynamical equation ( 6.22 ) in terms of $. Most of the terms in this equation combine into 
the usual GR-type equation for the potential <£. Few terms remain, however, and these are 
the contributions due to non-metricity. It is most convenient to write them in terms of X- 
We get: 

4TiGa 2 T5£ = - c 2 s A<S> + 371(1 + c 2 ) + [m' + H 2 (l + 3c 2 )] $ 



+ GrrGa 2 



5S" + m (5S' + H8S) - \a8S 



- ^f(P + P)Ax - ^G (p' - c 2 p') X ' ■ (6.25) 

Here the part which does not contain x 1S the- standard GR dynamical equation for the 
potential <]?, while the part containing x 1S the modification due to non-metricity. In 
obtaining ( |6.25D , we used the background field equations. Let us also give an expression 
for the constraint ( |6.18 ) in terms of the potential $. Again using the background equations, 
we obtain 

4TrGa 2 5p = A$ - 3U& - 37i 2 $ + QirGa 2 [ASS - 3U (SS' + HSS)) - AirGp'x • (6.26) 

The last term here is the contribution due to non-metricity. One can also give an equation 
for the velocity perturbation in terms of the new potential: 

AirGadu = + H<S> + QirGa 2 (5S' + 7i5S) - 4vrG(p + p)x ■ (6.27) 

Finally, we will also need an expression for equation ( |6,1S| ) in terms of the new potential 
$: 



ATrGa 2 5p = + 3H& + (W! + H 2 )<S> + QirGa 2 



5S" + m (ss' + nss) - ^ass 

-^-(p + p)A x -^Gp'x' ■ (6.: 
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6.6 Applicability limits of the theory 

As we have discussed above, our classical theory is applicable in the case of small deviations 
from the "metric" behavior, which is quantitatively expressed as ( |2.25| ). Since we have 



(2.24) and, in the linear theory, according to (5.6), 

1 f 1 



H tJ ss SB 13 



we have the following condition: 

k 2 , 



k 2 



< 1. 



.29) 



(6.30) 



a 2l/VI a 2 4nGK + 2n 2 /a 2 
This condition is not stronger than the usual condition of smallness of density contrast: 

2k 2 



5p + 3(U/a)5u-^AS 



P 



< 1 



3U 2 - A a 2 



1*1 < 1, 



(6.31) 



provided the denominator 4ttGk + 2H 2 in ( |6.30| ) does not approach zero sufficiently closely. 
We have to assume that this never happens and that always 



y2 

2hGk + —r > . 



(6.32) 



Note that this condition is also necessary in order that the change of variable ( 6.23j ) be non- 
degenerate. Recalling the definition (|4.4| ) of k, we obtain a constraint on the parameter 
(3 or g by considering a dust-dominated universe: f3 > —4/3, or g > —1. Considering 
the inflationary universe with p ~ —p, we get a stronger condition f3 > —1/3, or g > 0. 
Finally, considering the "stiff" matter withp = p (which is realized, e.g., by the kinetically- 
dominated regime of a scalar field), we have the constraint f3 < 2/3, or g < 1. Thus, the 
physical values of g or, respectively, /3 = g — 1/3 lie in the domain 



0<5< 1, 



1 -<^< 2 - 
3 H 3 



.33) 



In a similar way, we obtain a constraint on the parameter 7 by considering a universe 
which is dominated by the cosmological constant: 7 > —4/3. An even stronger condition 
can be obtained by considering a radiation-dominated universe. We then have: 



7 > 0. 



(6.34) 



We, therefore, assume that inequalities ( |6.33| ) and ( |6.34| ) are satisfied by a large margin. 
In particular, the quantity f3 cannot be much larger than unity by absolute value. 



7. Evolution of scalar perturbations 
7.1 Inflation 

The primordial spectrum of perturbations is obtained in the inflationary paradigm. There- 
fore, we start with application of our theory to the simplest model of inflation based on a 
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single inflaton field <p. We assume that the inflaton has the usual coupling to the metric 
defined by the £>-field. Then its action can be written as 



S\ 



(7.1) 



where g^ u is the metric defined by B. 

The background field equations for ip take the usual form 

1 12 , 2rrf J , A()a 2 

-(p +a V(p) 



8ttG 



+ 



U 2 -%' = AirGp 



r2 



(7.2) 



where, as before, the prime denotes the derivative with respect to the conformal time. 
The equations for perturbations can be obtained from the general equations ( |6.18 ) and 



( 6.20 ) by substituting appropriate expressions for 5p and Su: 



Sp 



Sp 

V 



2% 



5(p 

V 



From these expressions, we find 

op H du 

a 



S<p 

V 



GR 



+ n 



5u 



5p . 



dp 

V 



GR 



(7.3) 



(7.4) 



To obtain a dynamical equation for perturbations, it remains to express this quantity in 
terms of the potential <I> using ( gjg ). For this purpose, we first use equation ( p. 13] ) to 
express <I?gr in terms of the potential ^gr: 



GR > 



(7.5) 



where we have used expression ( |7.3| ) for 5u as well as the background field equations (|7^). 
Substituting (|7j) into (Q and once again using the background field equations, we can 
write equation ( |6.24j ) as follows: 

n 2 -H' f n 8p T V H 2 -H\ , A ^ 

^ J 7 + ^gr - „_,„, A X ' = A$, (7.6) 



3a 2 % 



where, as before, $ is given by J6.23[) . 

Now we use expressions (|6.20|) and (|6.10|) , respectively, for Su and ^gr in terms of x 
to transform the expression in the brackets of ( |7.6|) into the following form: 



H^ + ^GR 

iff 



n 2 /a 2 



W -W \ H 



a 2 $V 1 A 

+ 3^ AX ' 



(7.7) 



where we have used the definition (6.23) of the potential <£. Combining equations 
and ( |7.7D , we see that the term containing Ax' cancels. The remaining term with Ax can 
be converted to a term with A$ by using the relation ( |5.23 ) between x an d 3>. The final 
equation for the potential is 



u 2 -w 


' H 2 /a 2 /a 2 $\'~ 




u 


U 2 -H'\H J 





1 + 



u 2 -w 



3 (H 2 + GirGa 2 K) 



A$. 



(7.8) 
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Let us also present it in terms of the derivatives with respect to the physical time: 



H 9t 



H \ ( a ® 



H 



3 (H 2 + GttGk) 



(7.9) 



Apart from the term on the right-hand side proportional to H , this is the usual equation 
for the evolution of the relativistic potential during inflation. Note that the denominator 
H 2 +6ttGk in the last term on the right-hand side of (7.9) is positive by virtue of constraints 
( 103]) . 

The quantity QttGk in the case under consideration is given by the expression 



GttGk = 2irG(3g - l)(p 



A = (3<? - 1) 



3H 2 + -H 



+ 39A ° 



(7.10) 



As we have already noted in Sec. |2.q , the equations for perturbations, in general, are not 
invariant under the simultaneous change 



A n 



Ar 



A, V(<p) + 



A 



8vrG ' 



(7.11) 



which would be the case in general relativity. However, they are invariant after an addi- 
tional change 

[ i 

(7.12) 



1 1 , X 



In particular, such a transformation will arise in the case where the parameter £ 2 is fixed 
to be I 2 = 1/Agh.Q. 

During the inflationary epoch, we have <C H 2 , and one can neglect the term 
containing H on the right-hand side of ( |7.9| ). Thus, we can argue that the evolution of 
inflation-generated perturbations is unchanged in the theory under consideration, and we 
can take the standard flat spectrum for the potential $ as initial conditions. This conclusion 
for the generated spectrum can be justified as follows. Equations ( |6.24 ) and (7.4) imply 
the relation 

4-KGip' 2 ( a5ip s 



A3> + 47tGV' 2 $ G r 



(7.13) 



On spatial scales smaller than the Hubble scales, with the comoving wavenumber k satis- 
fying k 2 3> Ti 2 , we can neglect the self- gravity of the scalar field and quantize the inflaton 
perturbation 5<p on the background of a homogeneous and isotropic inflationary universe, 
obtaining the standard spectrum for the modes <5</?k- Then we can use equation ( 7,13| ) to 
match the quantities $ and 5ip on small (sub-Hubble) spatial scales during inflation. In 
doing this, the term 47rG<^ /2 $GR i n ( 7.13[ ) can be neglected (just as it is the case with a 
similar term in GR) since, according to fl6.11| ), it is estimated as 



4tt(V 2 $gr = 4ttGV' 2 [$ + O {H' 1 ) $' + O (?T 2 ) + O {U~ 2 ) A$] 



(7.14) 



and is small in view of the condition k 2 3> 7~L 2 3> 4-irGip' 2 . Thus, the quantity <1> on small 
scales will acquire the standard amplitude with scale-invariant spectrum and then will 
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evolve according to equation ( |7.9| ), which, as we noted, differs negligibly from the standard 
inflationary equation for the relativistic potential. 

In the regime of very long wavelengths, we can neglect all terms with Laplacians. In 
this case, equation ( [7.13 ) implies the relation 



$GR « - 

Substituting this into ( |7.3| ), we get 



a5ip 



5(p\ %5ip 
ip 1 ) ip' 



dp « -m — — = T <p5<p. 

Equations ( |6.10| ) and ( 5.11 ) with Laplacians neglected give 
^GR ~ $ + , $GR 



1 „ T~L i 
$ " "TX + "^X • 



Then, according to ( |7.2j ) and (7Jj), we have 



(7.15) 



(7.16) 



(7.17) 



(7.18) 



-2$ 



(7.19) 



where, in the last approximation, we have neglected the two terms — 4nG(f x' compared 
to which is legitimate during inflation. Comparing ( 7,16| ) and ( 7. 18] ), we get the usual 
relation 

dp 
P 

valid on super-Hubble spatial scales during inflation. Note that, in the same inflationary 
approximation, we have <I?gr ~ ^GR ~ $ on super-Hubble spatial scales. 

7.2 Evolution with a generic equation of state 

In this subsection, we would like to demonstrate that, on asymptotically large spatial scales, 
the relativistic potential <3> defined in ( |6.23| ) behaves just as in general relativity, including 
the cases where the universe experiences slow or rapid transitions between epochs with 
different effective equations of state. To see this, consider equation ( 6.25[) in the case of 
vanishing shear S for a system described by a generic equation of state p(p), which can 
interpolate between different regimes. Note that this assumption implies vanishing of the 
entropy perturbation, 5£ = 0, as well as of the last term in equation ( p. 25 ), which equation 
then takes the form 

- c 2 s A$ + m (1 + c 2 ) $' + [W! + H 2 (l + 3c 2 )] $ - -— (p + P )A X = . (7.20) 
By the standard change of function (see |Tl], Sec. 7.3]) 



u exp 



'l + <*)Udri 



equation (7.2C) is transformed to 



u 



-off 



An 



0" 



-u 



cx(p + p) 



0, 



1/2 



U . 



(7.21) 



(7.22) 
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where 







c eS 



2L 

H 2 



-1/2 



1 f P + PA 
a \ p + p 



1/2 



d + 



2(p + p) 



(7.23) 



(7.24) 



3P(p - 3p) + 3 7 p A + 4p 

and we have used definition ( 6.23; ) and neglected the cosmological constant contribution in 
the denominator. 

The long- wave solution of ( [7.22 ) is obtained when the term with Laplacian is neglected: 
where C and j]q are constants of integration. In this regime, the quantity 



3 p + p \ % 



(7.26) 



is conserved. 

Equation ( |7.22| ) can be recast in the form 



i) 2 



u 



4r^A ( 



(7.27) 



which, in particular, shows that the quantities 



u ( p + p \ a? ^ 

— = a u oc — <P 



and 



A - 



x9 + Ip+pa( 9+ y 



1 



(7.28) 



(7.29) 



3 p + p v ' n~ m 2 

should remain continuous during a rapid transition between different equations of state. 
Since a and T~L are obviously continuous in this case, it follows that <3? must also be contin- 
uous. 

The relations derived in this subsection imply that the inflationary and post-infla- 
tionary evolution of the potential $ on large spatial scales will reproduce those of general 
relativity. The only modification as compared with GR occurs on small spatial scales, 



where oscillations in ( |7.22|) proceed with effective speed of sound ( 7.24 ), different from c 2 . 
During some periods in the cosmological evolution, this difference can be small and lead 
to neglig ible effect in (|7.22|) . In this case, characterized by 



2(p + p) 



3P(p - 3p) + 3 7/9A + 4p 



<c c; 



(7.30) 



we will say that perturbations evolve in the regime of general relativity (GR) on all scales. 
In the opposite case, when the difference between c 2 s and c 2 is essential, we will say that 
they evolve in the regime of modified gravity. 
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From flOOl) and ( |6lTl) , one can obtain relations between the relativistic potentials on 
small scales: 



^ G r = $ + ^2 \ l A x + nx' 



$gr = $ + -4 ( ?A X - x" + Wx' 



A' 2 



12vrGa 2 K + 6% 2 



$, k 2 ^n 2 , 



1 + 



(3c c 2 ff - 2) k 2 
12irGa 2 K + 6H 5 



In the regime of general relativity, where 2ttGk S> H /a , we have 



GR 



1 



3a 2 



GR 



1 + 2 4 



2\ 



2e2 



(7.31) 



k 2 ^>% 2 . (7.32) 



(7.33) 



In deriving the above relations, we have taken into account that, in the high-frequency 
regime under consideration, x" ~ c egAx, so that the term proportional to x" m ( |6.11| ) 
gives a contribution of the same order as A^. We see that, even in the GR regime, the 
usual GR relation $gr = ^GR is violated for physical wave numbers k/a > l/l, i.e., on 
scales smaller than the scale of deformation i. 



7.3 Radiation domination 

The evolution of scalar perturbations in a universe filled with fluid are most easily analyzed 
using the system of equations ( 6.25 ), fl6.26|) . For simplicity, we set the shear to zero and 
restrict ourselves to adiabatic perturbations, 5£ = 0. 

At the stage of radiation domination, we have a = aorj, so that H = 1/r], H' = —H 2 , 
and Tl" = TH? . The speed of sound at this stage is c 2 = 1/3. Equation (|6.25| ) then 
simplifies to 



<£" _|_ 

V 



c e 2 ff Ad> = 0, 



where 



c eff 



7PA + Ap 



(7.34) 



(7.35) 



37M + 4p 

and, in using equation ( |6.23 ), we have taken into account that p = 3p, hence, k = 7Pa = 
const. Equation ( 7.34|) has the form of the usual equation for a radiation-dominated uni- 
verse apart from the fact that the effective speed of sound became time-dependent. 

In a very early universe, p S> 7/>a, and the effective speed of sound ( |7.35| ) becomes 
equal to the speed of light. This is the modified-gravity regime. In the course of time, if 
the quantity 7/9 a becomes dominant over p, the effective speed of sound (|7.35|) turns to its 
standard value l/\/3- This is the GR regime. As usual, the limit of GR is obtained by 
sending 7 — > 00. 

From the constraint ( 6,26| ), we find the expression for energy-density perturbations 
during this stage: 



A-KGa 2 5p = A$ - m& 



37PA 



37PA + 4p 



m 2 <$> 



4(4p) 2 



(3 7/ 9 A + 4p) 2 



(7.36) 
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where we have used the fact that k = 7Pa = const during radiation domination and that 
p' = —4Mp. Again, as 7 — > 00, this equation takes the standard GR form. 

As usual, the qualitative behavior of solutions of ( 7.34| ) is easy to understand in the 
two asymptotic regimes: kr] <C 1, corresponding to modes with wavelength larger than 
the Hubble radius, and kr] S> 1, corresponding to modes that already entered the Hubble 
radius. In the first case, the term A$ can be neglected, and 3> = const is a non-decreasing 
solution. In this case, the last term in equation ( 7.36| ) is dominant and gives a 4 5p oc 
with coefficient of proportionality being a function of time. Note that the usual general- 
relativity relation 5p/p = — 2$ is no longer valid in our theory. Indeed, this will hold only 
in the GR regime, where p <C 7Pa 5 while in the early universe we may have p S> 7PA- 
Thus, for kr/ <C 1, we have 



~ J 6$, /j > 7 pa (modified-gravity regime) 
p I -2<3? , p <C 7 pa (GR regime) . 



(7.37) 



As a consequence, the long- wave gauge-invariant energy-density perturbation 5p/p passes 
through zero for 7pa ~ P, which is an interesting phenomenon connected with the gauge 
choice in this model. 

Note that the standard relation 



5p 
P 



-!- « -2$ 



GR 



(7.38) 



remains to be true in both regimes, in view of relation ( [7.17 ) between the potentials $ 
and $gr valid in the long-wave approximation. Since it is the quantity $ which remains 
constant on large scales during radiation domination, both 5p/p and ^gr change in the 
transition between the regimes of modified gravity and GR. Relations ( [7.17 ) in this case 
read 



GR 



$GR 



1 + 



2(4p)< 



(37pa + 4p) 5 



3$, p> 7 pa (modified-gravity regime) 
<3? , p <C 7Pa (GR regime) . 



1 + 



4(9 7 p A -4p) (4p) s 
(3 7 p A + 4p) 3 



-3$ , p 3> 7pa (modified-gravity regime) 
3> , p >C 7Pa (GR regime) . 



(7.39) 



(7.40) 



For modes with krj S> 1, the term A<3? in ( 7.34j ) is important. The dominant solution of 
Q7.34 ) oscillates with frequency determined by the effective speed of sound ( |7.35| ) and with 
a decaying amplitude oc r]~ 2 c ^ 2 oc a~ 2 c ^ 2 . In this case, the first term on the right-hand 



side of (7.36) is dominant, and the relation between 5p/p and the potential $ is the same as 
in GR. The energy-density perturbation 5p/p is now an oscillating function with frequency 



determined by ( [7.35 ) and with amplitude oc c e ^ 2 , which is constant both in the modified 



and GR epochs and growing during the transition. The relativistic potentials on such small 
scales are described by the general equations ( |7.3l|) -( f7,33|) . 

Thus, as in GR, perturbations do not grow during radiation domination. This is true 
both for the modified-gravity epoch in early universe and for the GR epoch at late times. 
However, the relation between the energy density fluctuation 5p/p and the main relativistic 
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potential $ on the super-Hubble scales is different at these epochs, as shown by ( |7.37| ). The 
effective speed of sound ( [7.35 ) , which determines oscillations on sub-Hubble scales, is also 



different at different epochs. During the transition occurring for 7/? a ~ p, if it is realized at 
the radiation-dominated stage, the long-wave gauge-invariant energy-density perturbation 
5p/p passes through zero. 

Perturbations at the radiation-dominated stage are studied in greater detail in Sec. || 
below. 

7.4 Matter domination 

At this stage, we have p = 5p = 0, a(rj) oc rj 2 , and equation ( |6.25 ) simplifies to 

+ -cLA$ = 0, (7.41) 
V 



c 2 eS = ■ ( 7 - 42 ) 

6k + 4p 



where 



In this case, k = f3p + jp\ is not constant [j3 and 7 are defined in (|4 . 5[) ] . On sufficiently 
large spatial scales, k <C H/c e g, we have $ « const while, on small scales, k S> T-L/c e s, 
the potential $ oscillates with decreasing amplitude oc i] 3c c s^ 2 ^ a ~ 3 ^ 2c e s^ 2 anc ^ with 
effective speed of sound given by ( 7.42| ). This behavior is quite different from the case 



of general relativity, where the corresponding relativistic potential remains constant at all 
spatial scales. 

To determine the behavior of the density contrast, consider the quantity 

5p + 3(n/a)5u 
d m := , (7.43) 



which is related to <3? through (6.24): 

4TTGa 2 p5 rn = A$. (7.44) 

We see that, on large spatial scales, k <C H/c e g, the density contrast grows like in GR: 
5 m oc a, while, on small spatial scales, k 3> T-L/c c g, it oscillates with a slightly increasing 
amplitude oc aT l l 2 c C Q 2 oc a 1 / 4 and frequency determined by the effective speed of sound 



in (|7T42|) 



It is useful to obtain a closed equation describing the evolution of 5 m . Using the 
conservation equations (|6.16| ) and ( |6.17 ) as well as our assumption p = 5p = 0, we get 



p5' m = 3p (V GR + H<S> GR ) + ^ASu - K 2 ) 5u . (7.45) 



Taking into account equation ( 6.13 ) and the background equation H 2 — %' = 4nGa 2 (p+p), 



we see that the terms proportional to Su cancel, so that 



S' m = i [ Ax' + —A6u. (7.46) 
cr ap 
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Taking the second derivative of 5 m , we obtain 

81 = \ (A X " - 2HA X ') + A$ GR - - A5u , (7.47) 

or ap 

where we have used the conservation equation ( |6.17| ) once again. The following identity is 
a direct consequence of (^) and ( |6.10| ) : 

A$ GR + -jA (x" - HxO = yA 2 X + I^ttGk + A X . (7.48) 



Combining (|7.46| )-( |7.4£| ) and using the relation ( 7.44 ) between 5 m and Ax, we have 



C + - ^2 A ^ = *xGa 2 p5 m . (7.49) 



Finally, using the relation Q6.23 ) between 3> and x as we U as equation (7.44), we obtain 



51 + H6' m - 2p A5 m = AirGa 2 p8 m . (7.50) 

It is illuminating to rewrite the above equation in terms of the physical time: 

S m + 2H5 m 2P - A6 ™ = 4TrGp5 m , (7.51) 
6k + 4/3 a z 



where if = d/a is the physical Hubble parameter. Equation ( 7.51 ) allows for a straightfor- 
ward passage to the GR limit: one just has to send k —> oo. For finite k, we get an extra 
term that leads to oscillations on sub-Hubble scales, an effect absent in general relativity. 
In view of the definition k = (3p + 7PA> we can write: 

Cls = /3 + 4/3 + m/p • (7 - 52) 

Thus, at the early epoch with p 3> ~jp\ (if this condition is realized during matter 
domination), we have oscillations with effective speed of sound c 2 ff = 2/(3/3 + 4). If 
7 > 1, then, as the energy density decreases with the expansion, one may enter into 
the regime p\ -C p <C 7 Pa in which the usual GR behavior is recovered in a broad range 
of spatial scales. Note that the physical spatial scale on which oscillations take place 
k/a > H/c c ffd ~ o~ 3 / 2 (7Pa/ / o) 1 / 2 ~ l/l remains constant during the matter domination 
era. 

7.5 Lambda domination 

At later epochs in its evolution, the universe is dominated by the cosmological constant. 
It is thus necessary to consider the evolution of matter perturbations in such a universe. 
The scale factor in a Lambda-dominated universe is given by 
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so that % = — 1/77 and %' = l/f? 2 = T~L 2 - Equation ( 6,25| ) in this approximation reads 



- + 4 $ - c 2 ff A$ = , (7.54) 

7] 7] 1 



with the effective speed of sound 



2 AnGp 2p 



Ceff = 6ttGk + = (3/3 + 4), + (3 7 + 4) PA ~ 3^T4 (p/pa) < 1 ' (7 " 55) 

where we have taken into account that H 2 /a 2 = Ao/3 + 8irGp/3, k = (3p + "fpA, and the 
A-domination condition p^ 3> p. Equation ( [7.54 ) in this case differs from the general- 



relativistic analog only by the term with c 2 ff A<I>. 

On sufficiently large scales, k <C T-L/c e g, the evolution proceeds as in the case of general 
relativity, with the dominant mode 

$ oc 77 oc a' 1 . (7.56) 



The evolution of the density contrast is given by equation (6.26): 



4irGa 2 5p = A$ - ^Gp'x = A$ - ^c 2 s U 2 $ . (7.57) 

One can see that the evolution is the same as in general relativity, Sp/p oc a<£ = const 
in the range of scales c e s7i <C k 1i/c e s- However, the lower boundary in this case 
(Cg ff <C 1) corresponds to the scales well beyond the Hubble radius, which are unobserved 
in our universe. Thus, we may conclude that the long-wave perturbations, with k <C Ti/c e g, 
evolve in the GR regime. 

On small scales, k » H/c e ff, equation (|7.54j ) describes oscillations with adiabatically 
decreasing amplitude oc o,~ s ^ 2 c e ^ 2 oc a -3 / 4 (in the universe with pressureless matter). In 
this case, for c 2 ^ <C 1, we have 

4ttGcl 2 5p ps A$ , (7.58) 

so that the density contrast 5p/p oscillates with amplitude oc ar x l 2 c c ^ 2 oc a 1 / 4 , similar to 
what happens in the regime of matter domination. 

The threshold comoving frequency which distinguishes between these two types of 
behavior is growing as k\ := H/c e R oc a 5 / 2 , and so the physical threshold scale a/k\ 
decreases from its value ~ i at the matter-dominated stage as a/k\ oc a -3 / 2 . 

8. Effects of the modification of gravity on the power spectrum 

In this section, we give a qualitative analysis of the effect of modification of gravity on the 
linear evolution of matter and radiation density perturbations in the adiabatic mode. 



First of all, as we have seen in Sec. 7T| , the spectrum on super-Hubble scales generated 
during inflation is practically the same as in the usual inflationary theory based on general 
relativity. The subsequent evolution of the spectrum depends on the epoch at which the 
condition p = ^yp\ is reached, i.e., it depends on the value of the parameter 7. 
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Dealing with the post-inflationary evolution, we consider a two-component universe 
filled by non-interacting radiation and dark matter (we ignore the small contribution from 
the baryons at the radiation-dominated epoch). Using the basic system of equations ( 6.16) ), 
(6.17) and Q6.23 )— ( 6.28| ), it is possible to derive a closed system of second-order differen- 
tial equations for the two convenient variables, namely, the potential and the entropy 
perturbation 

3 

s m := -5 r — 5 m , (8-1) 

where, in this subsection, we have made the notation 5 r = 5p r /p r and 5 m = 5p m /p m for 
radiation and dark matter, respectively. One of these two equations is just equation ( 6.25| ), 
and the combined system in Fourier space reads [equation flS.3j ) is derived in the appendix] 

+ m (1 + c 2 ) + [k 2 c 2 ca + TH! + H 2 (1 + 3c 2 )] $ = ATrGa 2 c 2 p m s m , (8.2) 

1 h 2 n h A 

3c| Sm + %Sm + ~4^ Sm ~ 16vrGaV ' ( } 

System (|3.2j), ( |S.3j ) differs from the corresponding system in general relativity only by the 
presence of c 2 s instead of c 2 in one place in (^), where c 2 and c 2 s are respectively given 
by 

2 1 ( , , 3 Pm\ 1 2 2, 2 (P+P) 



Once a solution to ( j3^ ) and (|8 . 3[) is found, the density contrasts can be determined 
using (|87T1) and ( |6.26| ): 

s * - ~w~ r I 3 **' + ^ + 3 « 2 > •] + ^ + 3 ^ ■ < 8 ' 5 ' 

9c 2 3H 

6m = ~ 16^V ^ + ^ + 3 ^ 2) $] + " 3c ' Sm ' (8 - 6) 

where x 1S expressed through $ by ( |6.23| ). These equations formally differ from their 
general-relativistic counterparts only by the terms containing x- 



At the radiation-dominated epoch, equations (8.2) and (|8.3| ) take the simple form 

•qH" + 47/$' + k 2 c 2 ca r, 2 S> = ^Tr/sm , (8.7) 
tfs" m + V s' m + ^Trfsm = ^V<D , (8.8) 
respectively. Here T is a dimensionless parameter defined by 

(8ttG \ ^^ 2 

—Pe q a 2 eq ) , T ?? «l, (8.9) 

and a eq and p eq are the scale factor and radiation density at the radiation-matter equality: 

Pm — Pr — Peq- 

If the parameter 7 is such that the equality p = *yp\ is reached during inflation, then 
the evolution of the spectrum at the radiation-dominated stage is the same as in general 
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relativity. This follows from equations ( |8.2| ) and (p.q), which, in this case, coincide with 
those of general relativity at the radiation-dominated stage. Since inflation typically ends 
at energy densities p- m R ~ 10~ 12 /47rG 2 , and since 



1 



7PA 



8ttG£ 2 ' 



7» 1, 



we have the estimate for £ in this case: 



< 10°£ P ~ 10" 27 cm, 



1. 10) 



.11) 



where £ P = G 1 ' 2 ~ 1(T 33 cm is the Planck length. Such small values of £ make the theory 
practically indistinguishable from the general relativity (which, we remember, is obtained 
in the limit £ — > 0). 

Assume that 7 3> 1 is such that the condition p = ^p\ takes place at the radiation- 
dominated epoch. Let us determine the range of values of the parameter £ for which this 
is the case. The energy density p at the matter-radiation dominated stage is expressed as 



P 



8itG 



[1 + zY + n r (i + z) 



Mo 

SvrG r 



1 + 



l + z 



l + z, 



cq 



.12) 



where z eq is the redshift of matter-radiation equality, Hq is the current value of the Hubble 
parameter, and Q r and Vt m are the radiation and matter density parameters, respectively. 
When making numerical estimates, for definiteness, we will use the recently determined 
values of these parameters [O: 



n m /i? 5 «0.24, 



3 x 10 3 



.13) 



where /175 = //o/75km/s Mpc 1. 

The condition ^p\ = p defines the redshift zt : 



l + z e = 

and the condition zi > z cq implies 



l + z, 



cq 



3Hq Q m i 



1/4 



Z£ > Z, 



cq ; 



< [6n m H$(l + z eq ) 



31-1/2 



20kpc. 



5.14) 



.15) 



For the modes that enter the Hubble radius well in the GR regime, i.e., at z < Zg, the 
spectrum of perturbations will not be modified by the end of the radiation-dominated stage 
as compared to general relativistic cosmology, which can be seen as follows. Before the 
Hubble-radius crossing, the terms with k in system ( |8.2[ ) and ( |8.3[ ) can be dropped, after 
which this system coincides with the corresponding system in GR. After the transition to 
GR regime, when c 2 ff becomes equal to c 2 , the system again coincides with that of general 
relativity for all scales k, as was already discussed in Sec. [7J^. Therefore, the behavior of the 
solutions and s m is general-relativistic during the whole evolution. The non-standard 
terms with x m equations ( |8.5[) and ( |S.6| ) also become unimportant in the GR regime. 
Hence, the density contrasts are expressed through the basic functions $ and s m also in 
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a general-relativistic way. Therefore, their amplitudes and phases will not be modified 
compared to the general-relativistic expressions. 

The boundary value of the comoving spatial scale Xi = ao/kg that we were talking 
about is given by the condition that the corresponding wave crosses the Hubble radius at 
redshift Z£. We have 

h = Hl = TTi ' (8 " 16) 



from which, using ( 8,14| ), we get 

fa H V 1 + z. 



eq 



20kpc 



~ 100 



20kpc 



1/2 



Mpc, (8.17) 



where A eq = ao/k cq — 100 Mpc is the comoving spatial scale corresponding to Hubble- 
radius crossing at matter-radiation equality, and, we remember, this estimate works for 
£ < 20kpc. For A > A^, or k < = ao/A^, the spectrum of perturbations will not be 
modified by the end of the radiation-dominated stage. 

Consider now the modes that enter the Hubble radius at the modified-gravity epoch, 
where c^ s « 1. Such modes satisfy k S> kf, their comoving spatial scales, therefore, are 



considerably smaller than ( 8.17 ). In studying the evolution of adiabatic perturbations on 



such scales, for not very large values of rj, we can neglect the terms with T in equations 
Q8.7D and ( |8.8| ). Then, before the transition to the GR regime, we have c^ s ~ const = 1, and 
the exact solution of this system that describes the non-decaying mode with the property 
s m — > as rj — > is given by 



3$r 



sin kc e gr] 
(kc eS r]) 2 V kc eS r] 
$o |~cos kc e g7j — 1 



L eff 



-eff 



cos kc e gr] 
o 



.18) 



COS X 



-dx 



cos kc e grj — 1 



+ C + In kc e ffT] — Ci(kc e grj) 



3.19) 



where $ ~ $o = const at kr] <C 1, C ~ 0.577 is the Euler constant, and Ci(x) is the cosine 
integral. We have retained the constant c e ff in equations ( |8.18j ) and ( 8.19[ ) which, therefore, 



will be applicable also to the case of general relativity, where cj s 



1/3. 



Consider now a wave with a fixed comoving wave number k. It enters the Hubble 
radius at the moment r/ ~ Thus, when we are in the modified gravity regime 

the moment of Hubble-radius crossing occurs for a given wave earlier than would be the 
case in GR. The oscillations in $ thus start earlier, and the amplitude of the gravitational 
potential $ also drops more than in GR. A similar suppression effect is also present in 
the matter power spectrum. From the above expression for the entropy perturbation we 
learn that s m oc (kr]) 4 for kc e sr] <C 1. Thus, the entropy perturbation grows till the mode 
enters the Hubble radius where oscillations start. Since in the modified gravity epoch the 
Hubble-radius entry occurs earlier, the entropy perturbation grows considerably less then 
would be the case in GR and a suppression of matter power spectrum ensues. These effects 
can be estimated as follows. 
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Well after the Hubble-radius crossing (at kc e ff7] > 1), we can proceed to a new function 
u = t/ 2< 3?, for which equation flS.7|) with zero right-hand side gives 

u" + (k 2 (? eS - ^\ u = 0. (8.20) 

Therefore, an approximate solution to <3> in the regime k 2 c 2 s S> T~L 2 = r^ 2 can be given in 
the WKB form 

$ oc cos / kc e gdr] . (8-21) 

Ty/CeS J 



Comparing this with the leading terms in (8.18j), (|3.19[), we conclude that the WKB solution 



well after the Hubble-radius crossing can be approximated by 




cos 
c c ff Jo 



kctfffldr/, (8.22) 



^ cos f kcrt(r/)drf + (c - ~ + In kc° cfiV ) , (8.23) 

where c® s is the asymptotic initial value of the effective sound velocity, which we have 



2(<&Vc- 



inserted here for comparison with the case of general relativity. The second term in (|8.23| ) 



is just a source-free solution of equation (8.8); therefore, we have retained its form without 
any modification. 

We note that the same expressions (|8.22j ), (|8.23| ) will be obtained in general relativity 
with the substitution of c s for c e s everywhere. Therefore, after the transition from the 
modified-gravity regime, where c 2 s (rj) « 1, to the GR regime, where c 2 s (rj) = c 2 (rj) ~ 1/3, 



solution ( |8.22j ) for the potential $ will evolve just like in the case of general relativity 
except that its phase will be shifted and the amplitude of its oscillations will be lower by 
a factor of 

c -3/2 = 3 3/4 ^ 2 3 _ ^g_ 24 j 

This comprises the suppression by c 2 due to an earlier entry into the regime of acoustic 
oscillations with decaying amplitude, as well as amplification of the amplitude by the factor 

— 1/2 

c s that occurs in the transition from the modified-gravity to the GR regime. 

The amplitudes of the radiation and matter energy densities on small scales (kc e g S> H) 
will be given, according to (|8.5|) , (|8.6|) , by the expressions 




cos 



\-KGa 2 p r ( c ff ) 2 y C eff J Q 



kctf(rf)dr/, (8.25) 



~5 r - 3c s s m px - - - ^ 2 \l — [ 1 - ^3- ) cos / kc eS {r] )dr] 




(8.26) 
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We see that the amplitude of 5 r , compared with its general-relativistic counterpart, 
after transition to GR regime (where Cg ff = c 2 = 1/3), will be lower by the same factor 



(8.24), which will affect the CMB pattern on small angular scales. If one wishes to exclude 



this modification in the CMB spectrum at multipoles I < 2500, corresponding to the 



values reached by ACBAR [13], then one needs to impose the following constraint on the 
parameter I: 

l rec X e Xg A eq , / £ \ 1/2 200 „ . , , nnm . 

— -0.4 — < — — => ^<lkpc. (8.27) 



h A rec A eq A rec V 20 k P c / ~ 2500 

Here, Z rec ~ 200 is the characteristic multipole number corresponding to the Hubble radius 
crossing at recombination, and A rec = ao/k rec ps 250 Mpc is the corresponding comoving 
length. 

The amplitude of matter perturbation on small scales under consideration is given by 



.26 ). In general relativity, we would always have c^ s = c 2 s = 1/3, and the oscillatory 
term would always be absent. In the present case of modified gravity, the matter density 
contrast, apart from the usual monotonic evolution described by the last term in ( 8.261 ), 



also exhibits oscillations. However, these oscillations cease after transition to the GR 
regime (where c^ s = (? s = 1/3), and the matter density contrast is then described by the 
usual monotonic term. Its amplitude is smaller than that of general-relativistic expression 
by the factor c 4 = 1/9, although the argument in the logarithm is larger by a factor of 
v3- We thus conclude that, in the eventual power spectrum P(k) oc |<5 m (£;)| 2 of dark 
matter, there will be an additional suppression factor on small scales generated during the 
radiation-dominated stage: 



-Pmod(^) 

-Pacdm(^) 



C — 5 + In krj, 



cq 



9(C-i+ln^) 



2 

ps — , k » k e > k cq , (8.28) 
oU 



where P mo d(k) and Pacdm(^) are the power spectra in the modified gravity theory under 
consideration and in general relativity, respectively, r] eq is the conformal time at the matter- 
radiation equality, and k# and its corresponding comoving length are given, respectively, 
by ( |8.16| ) and ( p. 17 ). Such a strong suppression will place additional constraint on the 



value of t. In order that it does not disturb the observed power above the comoving spatial 
scales of ~ 1 Mpc, one would require Xi < 1 Mpc, or 

^<2pc. (8.29) 

On the other hand, suppression of the power of linear perturbations on comoving scales 
below 1 Mpc may be interesting from the viewpoint of the missing-satellite problem (see 
II). 

During matter domination, there will be an additional modification of power caused 
by oscillations on scales k > k m = %/c c k ~ aj\pil with slightly increasing amplitude 
5 m oc a 1//4 (see Sec. f^4| ). The corresponding comoving length scale is 

A„, = ^^(l + 2 ).8(Jil)(^) kp ,, (8.30) 
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At present, the evolution of perturbations on such small spatial scales is non-linear, and 
the linear-theory analysis of this paper is no longer applicable. On the other hand, it is 
not easy to probe the spectrum on such scales in the linear regime in the early universe. It 
is thus hard to see whether the new effect of matter density oscillations could be detected 
with the currently available data. 

The analysis made in this section, in particular, suggests that, if the condition p = 
7/j A took place at the matter-dominated stage, which requires values of I larger than 
those given in ( |8.15| ), this would modify the power spectra of dark matter and radiation 
very significantly compared to the case of general relativity, which is important from the 
viewpoint of current observations. For instance, the pattern in the CMB power spectrum 
would be affected because of the essential difference in the effective speed of sound 
from the usual one c 2 s in this case, and the matter power spectrum would be suppressed 
on all scales because of specific oscillations on the matter-dominated stage. Whether it 
is possible to fit all available observations in this case by varying cosmological parameters 
requires additional investigation using numerical methods. 

Finally, we note that we have studied here only the adiabatic mode, for which s m — > 
asymptotically in the past. For the bound ( |8.29| ) on the parameter I that arose in this 
section, the entropy mode is practically excluded, just as it is the case in the ACDM 
model. It also remains to be seen whether one can loosen the restrictions on the parameter 
I in the presence of an entropy mode. 



9. Discussion 



The modified theory of gravity that was studied in this paper has several specific features 
which make it an interesting object of investigation. Expressed in terms of the self-dual two- 
forms of the Plebahski formalism, the theory is local and has the same number of degrees 
of freedom as general relativity. Due to this property, the theory, in particular, respects the 
analog of the Birkhoff theorem, possessing a unique one-parameter family of (modified !) 
spherically symmetric vacuum solutions [||], which also turn out to be static. Unlike some 
other modifications of general relativity, such as massive gravity, it has a continuous limit 
to general relativity and, therefore, may be regarded as a smooth deformation of the latter. 
At the linearized level, the modification is described by a parameter I with dimension 
of length and a dimensionless parameter g; they enter the fundamental potentials V and 
R defined in ( 2.1 0| ) and ( p. 12 ), respectively. The value of g is restricted to lie between 
zero and unity from the requirement of absence of singularity in the equations for linear 
perturbations, while the parameter I is essentially free. The limit of general relativity is 
obtained as I — > 0, irrespective of the value of g. 

The cosmological properties of the theory also reflect the absence of new degrees of 
freedom. Thus, the dynamics of an ideally homogeneous and isotropic universe described 
by the modified gravity under consideration is absolutely the same as in general relativity. 
For a realistic universe, which is homogeneous only statistically, there will be the problem 
of averaging which might lead to modification of the effective Friedmann equations, the 
issue that lies beyond the scope of the present paper. 



- 36 - 



The theory of linear perturbations is modified in a very interesting way. Concern- 
ing perturbations of the scalar type, first of all, there exists an analog of the relativistic 
potential in this theory, which is related to the matter perturbations in the usual way 
through Eq. ( |6.24 ). Just as in general relativity, it is also possible to obtain a system of 



second-order partial differential equations describing the evolution of this potential; see 
Eqs. ( |6.23| )-( 6.28| ). This system is different from its counterpart in general relativity. Its 



analysis, performed in Sees. [?] and |8|, reveals the following basic qualitative features: 

1. The inflationary theory of generation of primordial perturbation remains intact be- 
cause corrections due to modification of gravity at the inflationary stage are sup- 
pressed by the small ratio \H\/H 2 <C 1. 



2. The evolution of perturbations at a radiation-dominated stage, described by ( 7.34| ) 



proceeds with a modified effective time-dependent speed of sound c e g , given by (|7.35| ). 
The relation between the energy density and relativistic potential, Eq. ( [7.36 ), is also 



modified as compared to general relativity. For sufficiently low values of the length 
parameter, i < 20kpc, given by ( |S,15| ), the transition from the modified-gravity 
regime (where c e g ~ 1) to the general-relativistic regime (where c e g ~ c s ~ l/y/3) 
takes place at the radiation-dominated stage. The modes that enter the Hubble radius 
after this transition evolve just as in general relativity, so that the spectrum in this 
long-wave region is not modified. The modes that enter the Hubble radius before this 
transition experience some suppression in the course of the transition; their amplitude 
drops by approximately 2.3, see Eq. ( |3.24 ), This will affect the CMB angular power 



spectrum on sufficiently small scales; in order that this spectrum remain unmodified 
below the multipole number I ~ 2500 reached by ACBAR |jl3| , the fundamental 
length parameter of the theory is further restricted to be I < lkpc; see Eq. ( |3.27| ). 

Modification of gravity during the radiation-dominated stage also has interesting ef- 
fect on the evolution of dark-matter perturbations on small scales, which is described 
in Sec. M. Specifically, the matter density contrast, apart from the usual monotonic 



evolution also exhibits oscillations, described by the first term in (8.26), that cease 
after transition to the GR regime (where c^ s = (? s = 1/3). Moreover, the amplitude 
of the monotonically growing part of matter perturbation is smaller than the similar 
expression in general relativity by the factor = 1/9. This leads to a suppression 
of the eventual power spectrum of dark matter by a factor ~ 1/80 [see ( |8.28|) 1 on the 
scales that enter the Hubble radius well before the transition to general-relativistic 
regime. In order that this suppression does not disturb the observed power above 
the comoving spatial scales of ~ IMpc, one would require I < 2pc, which turns 
out to be the strongest constraint on this parameter in the present scenario. On 
the other hand, suppression of the power of linear perturbations on comoving scales 
below 1 Mpc may be interesting from the viewpoint of the missing-satellite problem 
(see 0). 

3. The evolution at the matter-dominated stage is qualitatively different from that of 
general relativity because equation ( [7.41| ) for the relativistic potential still contains a 
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nonzero effective speed of sound, given by (|7g|) or ( ff752|) . For this reason, the closed 



equation ( [7.51 ) for density contrast exhibits oscillatory character on sufficiently small 



scales; the corresponding comoving spatial scale is given by (|8T3C| ). In view of the 



previous constraint £ < 2pc, this behavior will only affect the comoving spatial 
scales smaller than about 8kpc. Note, however, that, at present, the evolution of 
perturbations on such small spatial scales is non-linear, and the linear-theory analysis 
of this paper is no longer applicable. On the other hand, it is not easy to probe the 
spectrum on such scales in the linear regime in the early universe. It is thus hard to 
see whether the new effect of matter density oscillations could be detected with the 
currently available data. 

We should note that the above description is of qualitative character and, therefore, 
needs to be substantiated by numerical integration of the derived system of exact equations 
for perturbations, which will be the subject of subsequent work. We also leave to the future 
work an analysis of the effect of modification on the gravitational waves, with a preliminary 
description given in Appendix [B|. 
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Appendix 

A. Derivation of equations describing the radiation— dark-matter system 



The conservation equations ( |6.16| ) and (6.17), written in Fourier space for a system of 



matter and radiation with vanishing shear, read 



5> m = e m + 3V GR , 6' r = ^e r + w' GR , (A.i) 
-ue m - k 2 $ GR , e' r = -k 2 <s> GK - ~k 2 s r . (A.2) 



' rn 



Here, the variable 9 a for each component is defined as 

. k 2 5u a 



Pa + Pa 



(A.3) 



Equations ( |6.26 ) and ( |6.27| ) in the case of vanishing shear can be expressed, respec- 
tively, as 

4vrGa 2 (p m S m + p r S r ) = -2>U ($' + H§) - k 2 § - AttGp'x' , (A.4) 
AnGa 2 (p m m + ^p r e}j = —k 2 [$' + U<S> - 4ttG(p + p) X '] ■ (A.5) 
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Differentiating the entropy perturbation s m = |<5 r — S m , using (|A.1| ) and ( A.2), we get 



-8 r . 



(A.6) 



It remains to calculate the left-hand side of (|3lj|). Using the definition (JO|) of c 2 , we 



have 



—s m + Us m + -—s r , 



1 + 



3p w 
4/9r 



A: 2 



& + ft 



+ 



4 p r 



-S r — S r , 



~. ( PrrSm + ^Pr^r ) ~ ~. (Pm$m + Pr^'. 



4p r 



k 2 



m k 2 



4p r 



[& + n$-47rG(p + p)x] 



+ 



4p r 4irGa 2 



(A.7) 



Here, in the third equality, we have used equations flA.4j ) and ( |A.5| ), and, in the last equality, 
the conservation equation p' = —3'H(p + p). Equation ([A .7]) is the desired result (|S,3D. 



B. Tensor sector 
B.l Curvatures 

To describe the linearized connections and the curvature in the tensor sector it is very 
convenient to introduce the following two differential operators on the space of symmetric 
traceless transverse matrices: 

DX ij ._ X tk^jkl + ^ X ijy ^ Dx ij ._ X ik i£ jkl _ ^ X ijy ^ (B.l) 

where, as usual, the prime denotes the conformal time derivative. It is easy to see that the 
operators D and D transform symmetric transverse traceless matrices into matrices with 
similar properties and that 

DD = DD = -□, (B.2) 

where we have denoted □ := — d 2 + A. In terms of the derivative operators introduced, we 
get 

SAi = -i,D(ax ij )dx j , 5Ai = --D(ap^)dx j . (B.3) 
x a d y a 

For a linearized connection of the form 5 A 1 = A l] dx 3 with a symmetric traceless matrix 
A lJ , the linearized curvature is given by the following expression: 
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Thus, we have 



V Q 5A\ 



2a 



-D 



Vo6A* p = --D 



D (a X ij ) 



D (aft*) 



D (a X ij ) 



^0 ' 



(B.5) 
(B.6) 



Note that the derivative operators D and D do not act on the basis two-forms in these 
expressions. 

B.2 Field equations 

For tensor perturbations, there is no source on the right-hand side of Q4.S| ), and we get the 
following simple system of equations: 



aD 
1 



-D 



a 



D (a X ij ) 



D [a X ij ) 



1 



— {n 2 + %') X ij + aD 



D (apV 



{U 2 - W) p ij = 8TrGK X ij , (B.7) 



{U 2 - W) X tj + -DD(ap») - (n 2 + W) p lJ = . 



(B.8) 



Now expanding the derivative operators one obtains the following coupled system of second- 
order differential equations: 

(d 2 - A - 2Ud r] - m 2 - 8irGa 2 K + 4iHe) X ij = (d 2 + A + 2id v e) p ij , (B.9) 
a 2 (d 2 — A + 2Hd n ) p lj = (8 2 + A - 2id v e) X ij , (B.10) 

where e is a first-order differential operator that acts on the space of symmetric trace-free 
matrices via: 

eX ij = X ik ie jkl _ 

We note that e 2 = —A. As a check, in the limit x^ ~^ such that kx^ is kept finite, we 
get the GR result (d 2 - A + 2Ud v )p lj = 0. 

B.3 Action principle and Hamiltonian analysis 

We have written the system of equations ( [B.9D , ( B.10D in the above form to make it obvious 
that they can be obtained as Euler-Lagrange equations for the following Lagrangian: 

a 2 

Ctens = --P {d 2 -A + 2Hd v ) p + X {d 2 + A + 2id v e) p 



—^X (d 2 ~ A - 2Ud n - m 2 - SitGna 2 + 4i^e) X , (B.12) 



2a 2 

where we have dropped the internal indices for brevity. Indeed, one easily obtains both 
equations using the fact that the operators d 2 — A + 27id v , d 2 , and d 2 — A — 2718^ are 
self-adjoint with respect to the inner products given by 



dr/d 3 x a 2 <frip , (cf>, ip)o = I drjd i x (f)tjj 



,i/>)- = J drjd 3 x^(j)ifj, (B.13) 
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respectively. Moreover, the operator e is self-adjoint with any of these inner products. 

It is clear from the form of this Lagrangian that it is degenerate and thus describes 
just one propagating field, not two. Indeed, the time derivatives of the fields enter in the 
combination (a 2 /2)(p' — x'/a 2 ) 2 , which makes it clear that the system is degenerate. There 
are many ways to unravel the dynamics of such a system, but the most powerful, if not 
always the easiest method is via its Hamiltonian analysis. Thus, let us rewrite it in the 
Hamiltonian form. To this end, we first rewrite the Lagrangian in the form that makes 
it obvious what the canonical momenta are. Integrating by parts in the terms containing 
second time derivatives, we write: 

a 2 

Ctens = y {p' 2 + pAp) - (ep - ip')(e X + tf) 

+ ^ W 2 + X&X - Wi X ex + (4^ 2 + 8ttGci 2 k) X 2 } . (B.14) 
Now the conjugate momenta are: 

k p = a 2 p' - x' + iex , ^ x = ~2 X* ~ p' ~ ie P • (B.15) 

cr 

Thus, the following primary constraint holds: 

fa := vr p + a 2 ir x + ie(a 2 p - x ) = . (B.16) 

The Hamiltonian is obtained as "Htens = K p p' +n xX ' — C tC ns+ufa, where u is a Lagrange 
multiplier and <f) is the primary constraint. The Hamiltonian then must be written in terms 
of the momenta. Only one momentum is needed for this, for which we choose ir p . We get: 

%tcns = ^ ~ ~2V e X - Y pAp ~ pAx + 2a 2 ( 4i ^ xex ~ 4 ^ 2 ^ 2 ~ 8kGci 2 kx 2 ) + u<j>\ . 

(B.17) 

As a check, one can verify that the Poisson brackets of this Hamiltonian with the phase 
space variables ir p , tt x , p, and X §i ve r i se precisely to the original equations (|B.9| ), ( B.1C| ) 
with the Lagrange multiplier u = X ' l& 2 as a consequence of x' = {-^tensj x} = a2u - 

Now the condition that the constraint (j>\ = is preserved in time gives the secondary 
constraint: 

2 : = 0; = 2a 2 ^(^ x + iep) + {H tcns , fa} = . (B.18) 

We find: 

4> 2 ■= a 2 Ap + ievr p - Uit p - YHex + 2H 2 x + 4irGa 2 K X = , (B.19) 

where we have used the relation a 2 (7i x + iep) = iex — n p that follows from fa =0. As a 
check, we note that (B.19) is precisely the difference of the two original equations (|B.9|), 



(B.10) 



Now one can see that the Poisson bracket of fa with fa is a constant, and so the 
constraints are second class. These constraints can be used to solve for tt x and x i n terms 
of the phase space variables ir p and p. After that, one should substitute the result into 



( B.17 ) to obtain an effective Hamiltonian with second class constraints solved for. In 



addition, one also has to compute the arising Dirac bracket. 
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To solve flBl9D , we introduce: 



m 



2 := 2H 2 +4irGa 2 K. 



(B.20) 



We now get for \\ 



e — m 



,2\-! r 2 



[a 2 Ap+(ie-n)Tr i 



Pi ' 



(B.21) 



where the inverse of (p-Le — m 2 ) should be interpreted by passing to the momentum space. 
There is no operator-ordering ambiguities in this formula as all the operators only involve 
spatial derivatives and thus commute. 

We now note that Hamiltonian (B.17) can be rewritten on the constraint surface ( |B.1G| ) 

as 

/eff 



H 



1 a 2 % 



(B.22) 



where % 1S to be substituted from (B.21). The action in the Hamiltonian form on the 
constraint surface now reads: 



? eff 
'tens 



drjd x 



I . X • A. I ■ 2 \ n i 

p 7T p + — iex 2 V^P + lfl 6 P) ~ ™ 



X 



eff 
tens 



(B.23) 



where we have used the constraint ( |B.16| ) to express tt x in terms of the other phase space 
variables. Now, the term 



x' . H . f x^x\ 
- 2 *x - - 2 x*x = [^r ) 

here is a total time derivative and can be dropped. This leaves us with 

,2 



J, 



tens 



/ 



drjd x 



P >7T p -^2 ( < K P + ia2e P) ~ ^2 + Y pAp + ^ X7Tp 



(B.24) 



(B.25) 



where ( B.21 ) must be substituted. After this is done, and all the dust settles, we get the 
following effective action in the Hamiltonian form: 



eff 



drjd?x ( p 



1 



+ Ap^ 



m 2 -\He 2a 2 p 
a 



1 + 



U-ie 
m 2 — YHe 



\ 2 -U I +^7P A P 

a \ m z — lHe J 2 



i + 4 



a 2 \ m 2 — lHe 



(B.26) 



As before, the limit to GR is easily obtained by sending m 2 — > oo. Another simple limit 
is that of passing to the Minkowski spacetime background. This is obtained by setting 
T~L = and a = 1 everywhere, as well as taking m 2 to be time-independent. One again 
obtains an unmodified system, even for a finite m 2 , which is consistent with the earlier 
observations in the literature that the gravitational waves in the Minkowski spacetime are 
unmodified. However, we see that tensor perturbations around an expanding universe are 
modified quite non-trivially. 
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This action can be further rewritten in the usual form depending on generalized ve- 
locities only by integrating out the momentum ir p . We get 



a. 



eff 
tens 



drj (fix ■ 



1 + 



pAp 



a 1 



77V 



iUe 




-i 



777 



777' 



+ 



Ap 



m- 



lHe 



(B.27) 



The resulting effective action is complex. The prescription for dealing with such complex 
actions advocated in || was to require the metric to be real and take the real part of 
the action. However, a different prescription may also be possible. 4 Thus, it may be that 
the appropriate prescription is to allow the action to be a complex, holomorphic function 
of the metric g pu , with the latter treated as a collection of complex variables. Then, the 
condition that the imaginary part of the action vanishes gives a relation allowing to express 
the imaginary part of the metric in terms of its real part. The real part of the action can 
then be written as a functional of the real part of the metric. More work is needed to decide 
which prescription for dealing with "reality conditions" for this theory is appropriate. We 
leave an analysis of the tensor sector that depends on this prescription to future research. 

An alternative quick way to obtain the field equation for p that stems from ( |B.27| ) is 
as follows. One notes that equation ( |B.10| ) can be written as 



7r„ 



a 2 Ap 



ie(ie - d v )x = . 



(B.28) 



On the other hand, from the definition ( B.15[) of the momentum 7r p , it follows that 

(ie - d v ) X = vr p - a 2 p' . (B.29) 
This gives an equation that involves only p and tt p : 



a 2 p 



ie(7r p 



2 l\ 

a p) 



0. 



One can rewrite this as 



(d v - ie)(vr p + ia 2 ep) - 2Ha 2 iep = . 



(B.30) 



(B.31) 



Using the expression for 7r p that was obtained above (when integrating it out), or substi- 
tuting x into the definition ( |B.15| ) of the momentum tt p and solving for tt p , we obtain the 
following compact expression for the combination 



tt p + ia ep 



1 + 



n 



ie 



(d v + ie) 



m 2 p 



m 2 — iHe 

As a result, we get the following differential equation for p: 

U 



YHe 



(0, 



ie )a 



1 + 



ie 



777 z 



me 



-1 



(d v + ie) 



m 2 p 



2Ua 2 \ep = . 



(B.32) 



(B.33) 



777^ — me 

As a check, we note that the resulting equation takes the general-relativistic form (a 2 p')' — 

a 2 Ap = in the limit m 2 — > oo. We shall not analyze the arising modified gravitational 

wave equations in this paper, leaving this to future work. 

4 One of the authors (KK) is grateful to Laurent Freidel for a discussion that led to an alternative 
prescription. 
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